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Abstract 

We calculate supersymmetric QCD corrections (squark/gluino loops) to quark pair 
production in e^^e" annihilation, allowing for mixing between left- and right-handed 
squarks and taking into account the effects of nonzero quark masses. Corrections to 
the Z boson partial widths are generally small and positive, except in the case of large 
h squark mixing, where they become negative. At high-energy e^e~ colliders, larger 
corrections to the total cross sections are possible. Corrections to forward-backward 
asymmetries are negligible except possibly for top quarks, where they are sensitive to 
t squark mixing. We also comment on the possibility that the gluino mass is only a 
few GeV. 
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1. Introduction 



The introduction of Supersymmetry (SUSY) is one of the most attractive extensions of 
the Standard Model (SM). It not only stabilizes the huge hierarchy between the weak 
scale and the Grand Unification or Planck scale against radiative corrections; if SUSY is 
broken at a sufficiently large scale, as is the case, e.g., in Supergravity (SUGRA) models 
[Q, it might allow to understand the origin of the hierarchy in terms of radiative gauge 
symmetry breaking [Q. Moreover, SUSY models offer a natural solution of the cosmological 
Dark Matter problem |p, and allow for a consistent Grand Unification of all known gauge 
couplings, in contrast to the nonsupersymmetric SM [^]. All these attractive features are 
already present in the minimal supersymmetric extension of the SM, the MSSM, to which 
we will stick in this article. 



Unfortunately no direct signal for the production of superparticles has yet been observed; 
experimental searches so far have only resulted in lower bounds on sparticle masses, the 
most stringent ones coming from LEP [0] and the Tevatron ||^. It is therefore tempting to 
look for SUSY through precision measurements, where quantum corrections involving su- 
perparticles might alter SM predictions. The potentially largest corrections can be expected 
from corrections involving strong interactions, i.e. from squark and gluino loops. Given 
the inherent uncertainties of cross section calculations as well as measurements at hadron 
colliders, the most promising (and also the simplest) process where such corrections can be 
probed is quark pair production in e^e~ annihilation^ 

In this paper, we calculate the supersymmetric QCD corrections to quark pair production in 
e^e~ annihilation, allowing for mixing between left- and right-handed squarks and taking 
into account the effects of nonzero quark masses. At LEPl energies, we find that these 
corrections are small and positive for Z decays into light quarks; however for hh final states, 
mixing in the h squark sector can affect the correction to the cross section, and can even 
change its sign. In the case of top quark pair production at high-energy e'^e" colliders, 
the effect of mixing in the t squark sector on the total cross section is less significant, since 
the dominant photon exchange contribution is not sensitive to it. The correction to the top 
forward-backward asymmetry does depend on the details of t squark mixing but unfortu- 
nately the correction is always very small, and will therefore be difficult to measure. 



Supersymmetric QCD corrections to quark pair production in e^e~ annihilation were first 
discussed in Ref. for LEPl energies in the approximation of negligible quark masses and 
squark mixing and of equal masses of the superpartners of left- and right-handed quarks. 
In Ref. |Tl[] the effect of squark mixing has been included at LEPl energies and found to be 



small. However, in that paper only corrections to the Z-quark couplings present in the SM 
at tree level are considered, while we compute all CP conserving form factors for both the Z 
boson and the photon (the latter are needed for cm. energies away from the Z resonance). 

^Squark and gluino loops also contribute to rare K and B meson decays and oscillations. However, these 
corrections always involve flavour-changing couplings, which in the MSSM are induced only through weak 
interactions. As a result, in the MSSM supersymmetric QCD loops in K and B meson physics are actually 
smaller jo) than loops involving electroweak gauginos or Higgs bosons. 



In the limit of zero quark mass and squark mixing, our results for the total cross section 
fully agree with Ref.[|TU| both numerically and analytically; we also find general numerical 

Finally, we also compute corrections to the forward-backward 
focussed on corrections to total rates. 



agreement with Ref. [11 



asymmetry, while the previous papers ||T0|, |TT 



The rest of this paper is organized as follows. In sec. 2 we set up the formalism and present 
our analytical results for the corrections to the most general set of CP conserving '-yqq and 
Zqq couplings. In sec. 3 we show numerical examples both for LEPl and for a future high- 
energy e+e~ linear collider operating at y/s = 500 GeV. Sec. 4 contains a summary and 
some conclusions. For the convenience of the reader explicit expressions for the scalar 2- 
and 3-point functions appearing in our results are listed in the Appendix. 



2. Formalism 

The most general Zqq and 'yqq vertices compatible with CP invariance can be written as 



'X. 



where eo is the electric charge of the proton, P = pi — P2 with pi,p2 the momenta of the 
quark and anti-quark and g'^ = 1, g^ = I / (iswcw) with = 1 — = siia^ 9w- Because 
CP is conserved by strong interactions, terms proportional to Pfj.'j^ should be absent and 
this fact provides a good check of the calculation. In principle one can also have scalar and 
pseudoscalar couplings, and 5^75 where q = pi + p2 is the momentum of the gauge boson; 
but in e~^e~ collisions these terms give contributions which are proportional to the electron 
mass and are therefore totally negligible. At the tree level, S^''^ vanish, while the vector and 
axial-vector couplings take the usual form: 

(y/r = = 2/f - 44e, , (Af )° ^ a, = 2/f , {V^'f = e, , = 0, (2) 

with = ±1/2 the weak isospin and Cq the electric charge of the quark. When loop 
corrections are included, S^''^ terms appear and the bare vector and axial-vector couplings 
are shifted by an amount 

§Vf''' = Vf'^-iVg^'^)° , ^ = ^ - (Af ^)°. (3) 

In previous work [1T0| , p!T[] only the corrections to and were considered. We find 
that even for heavy (top) quarks the corrections coming from the scalar form factors S^''^ 
are indeed somewhat less important than the corrections to the couplings that are already 
present at tree level. 

Since we are interested in radiative corrections involving strong interactions, we only need 
to consider diagrams involving squark and gluino loops. As stated in the Introduction, we 



will include effects proportional to the mass of the produced quarks. As well known |12|, the 
supersymmetric partners of left- and right-handed massive quarks mix; the mass eigenstates 
qi and q2 being related to the current eigenstates and (Jr by 

qi = qL cos 6 + qRsin 6 , q2 = -qLsmO + qncosO. (4) 



The mixing angle 6 as well as the masses m^^, rriq^ of the physical squarks can be calculated 
from the following mass matrice^: 



Ml 



2 

rriz 



, mf + 0.35Dz -mt{At + ficotp) 
-mt{At + cot (3) m] +m^^ + Q.lQDz 



mf^ +ml- 0A2Dz -mf,{Ab + /itan f3) 
—mb{Ab + /itan/3) 



m? +mf-0.08Dz 



(5a) 
(5b) 



where Dz = Mz cos2/3, tan /3 being the ratio of the vacuum expectation values of the two 
neutral Higgs fields of the MSSM 'fnii^fn'bn breaking masses, A^^t are parame- 

ters describing the strength of nonsupersymmetric trilinear scalar interactions, and /i is the 
super symmetric Higgs(ino) mass, which also enters trilinear scalar vertices. Notice that the 
off-diagonal elements of these squark mass matrices are proportional to the quark mass. 
In the case of the supersymmetric partners of the light quarks mixing between the current 
eigenstates can therefore be neglected. However, mixing between t squarks can be sizable 
and allows one of the mass eigenstates to be much lighter than the top quark. Sbottom 
mixing can also be significant if tan/5 ^ 1; even in supergravity models with radiative sym- 
metry breaking tan/? can be as large as mt/mi, |jl3[. 



The interactions of the photon and the Z boson with squark current eigenstates are described 
by the following lagrangian [Q: 



qqV 



i=L,R ^wcw 



(6) 



After the introduction of nontrivial squark mixing, this becomes |]n 
q*i df,qi + q; d^q2 ^—Z^" 



^qqV — icA^Cq 



SwCw 



sin e cose{q{ 9^ + ql qi) 



+ (/f cos^ d - slreq)ql qi + {if sin^ 6 - slreq)q; 9^ ga) 



(7) 



Finally, the squark-quark-gluino interaction lagrangian in the presence of squark mixing is 
given by 



jCgqq = —iV2gsT°'q (cos9qi — sm9q2^ \ ~ (sin6'gi + cos 



75 



-g'^ + h.c, (8) 



where gs is the strong coupling constant and are SU(3)c generators. Note that in eq. 



we have assumed Mg > 0. 



Including the corrections due to the squark/gluino vertex diagram shown in Fig. la, and 
taking into account the mixing between the left and right-handed squarks as well as the 
finite mass of the external quarks, the photon couplings to quarks are shifted by: 



^We ignore generation mixing between squarks, which in case of the MSSM is only induced radiatively 
by weak interactions. 



3 TT 2 



Cl' + cl 



(5'j')vE 



3 TT 2 

4 a., e„ 



rill ^22 



+ C2 ) - ^gMg sin 20(^1^ - 



-Y22^ 



^'^^ 3 TT 2 

The corresponding shifts of the Z boson couphngs to quarks are: 
4a, 

3 TT 



(2/f cos^ ^ - 2sl,e,)Cl^ + (2/f sin^ ~d - 24e,)C3' 



4 a. 



^VE — 77 



3 TT 



(2/3^ cos^ ^ - 2s^e,) COS 2901^ - (2/^^ sin^ 9 - 2s^eq) cos 2^C|2 



+ /f sin2 2e~(Cf .^Cg^ 



(9a) 
(9b) 
(9c) 

(10a) 
(10b) 



('S'5)vE 



4a, 

'3 TT 



(2/3^ cos^ ^ - 2s^eg)(mJ(:7ii - m^M-g siyi29C\^) + (2/^^ sin^ 9 - 2s^eg) 

X {mlCf + m^Mg sin 2^Cf ) + if sin 2^ cos 29mqMg{Cl^ + )] . (10c) 

For 2, j = 1, 2 the C|-^2,3 cire defined as C^-' = Cfc(s, rriq, m^., mg^., M^). We use the Passarino- 
Veltman reduction to scalar integrals |]15|, and the intermediate function and C!f defined 
as (note that we use a slightly different notation than in |p!6|): 



-1 

2^ 



25o(s, mq^.TJiq^) - Boiml, Mg, nig^) - 5o(mq, Mg, m^-J 



+ (2M/ + 2mJ - _ 



= — 
" 2s 



Bo{ml, Mg, rrigj - Bo{ml, Mg, mq^) + (m^/ - mq^^)Cl^ 
In terms of these functions, the C]^ are given by: 



^i — ~ '^^+ 1 

4 



C^s' + ^K-.' + m,/ - 2M/ - 2ml)Ci - \bo{s, rrig^^m,;, 



+ \Bi{ml Mg, m,-J + -^B^{ml Mg, m^;, 



4 



2M/C^^' + 1 + 5o(s, m,., m,-J + (m,-,^ + m,-/ - 2m^ - 2M/)C: 



+ (m,,2-m,/)C!.^- 
Here, Bi is given by 

1 



Bi{s,mi,m2) 



s + m^- m2)Bo{s, mi, 7712) + Ao{m2) - v4o(mi) 



(11) 



(12a) 



(12b) 



(12c) 



(13) 



and the functions A^, Bq and Cq correspond to the scalar one, two and three point functions 
P!7| , respectively, and are given in Appendix. 



The renormalized vertices are derived by adding the counterterm originating from the on- 



shell self-energies of the external quarks. Fig. lb. Following the procedure outlined in |T^, |18 
one obtains 



CT 



{^Vg)c'T = Vq6Zv + a^SZA 



qjCT — GqSZA, 



{6A^)cT = aqdZv + VqSZA- 



(14) 



Here, 5Zy and 5Za are given by 
1 as 



SZ. 



V 



3 TT 



Si(m2, M^, m,-J + 5i(m2, Mg, m,-J + 2mJ((S;(m2 M^, m,-J + S;(m2 M^, mg-,; 



-2mqM~g sm2e{B'o{ml, Mg, m,-J - ^^(m^, M~g, m,-J) 



SZa 



1 

3 vr 



cos 26 



Bi {m^q , Mg , rriq^ ) - 5i (mj , , rriq^ ] 



(15 
(15 



The full SUSY-QCD correction to the vectorial and axial couplings is just the sum of the 
unrenormalized vertex correction and the quark self-energy counterterm^ 



= (5V^)vE + (5n''^)cT, 
5Af = (M^'^)vE + (5^^'^)cT. 



(16a) 
(16b) 



The expressions (^ - (|T5D are rather cumbersome. For many applications squark mixing 
can be neglected. If one in addition assumes approximate degeneracy for the squarks, rriq^ = 



m. 



Q2 



Trig, the corrections simplify considerably, and one finds: 



^ on 



5Ar = \-{ArfC, 
^ 6 n 



3 TT 



where C is given by: 



ml{V^'^fC2{s,mq,mq,mg,Mg) 



C = C^{s,mg,mq,mq,Mg) - ^Bi{ml,mg, Mg) - ^mlB[{ml,mg, Mg). (18) 

Furthermore, for massless final state quarks the Sg term vanishes and the correction to the 
axial and vector couplings can be expressed by a single two-dimensional integral as: 



(17) 



1 



2 Jo 



xdx / dy log ■ 



Ximg"^ -Mg'^)+ Mg'^ 



-sx'^yil -y)+ x(mg2 - M^^) + M^^ ' 



(19) 



•^Note that in this convention the vector couphngs Vq and scalar couphngs Sg do not reduce to their SM 
values even in the limit of infinite squark masses; of course, very heavy squarks and gluinos do decouple 
from physical observables such as cross sections. Physically equivalent results can be obtained by ignoring 
the diagrams of Fig. lb, and performing the renormalization by simply substracting the corrections at zero 
momentum transfer, s = 0; in this scheme separate counterterms for Vq and Sg can be defined, so that each 
coupling by itself reduces to its SM value in the limit of large sparticle masses. 



in agreement with |T^. For large squark and gluino masses, mg,Mg ^ s, the correction is 
just 



C 



1 



12 K-2 - M/)4 



(20) 



If in addition the gluino mass can be neglected compared to the squark mass, one simply 
obtains C ~ s/(36mg^). 

In terms of the vertices (|I]), the differential cross section da{e^e~ — ^ qq)/d cos 9 reads (we 
define 6 as the angle between the quark and the incoming positron): 



da 



dcos^ 8 



|iV,/5, {D,,el [(2 - /5J sin^ eWj f + f3'^{l + cos' ^)(A 



^^2_2/5jsin2^V7^g^ 



+ Dz,e,Ve [(2 - W ^)\/;\// + /5j(l + cos^ e)A^Af - /5j sin^ ^(\/;5f + 
+ Dzzivl + al) [(2 - (31 sin' d){Vg^Y + /3j(l + cos^ ^^)(Af )2 - 2(3^ sin' ^V/Sf 
+ 2Dz7eeae/?, cos eiVg^A^ + Vg^AD + %DzzaeVel3, cos eVg^A^^} . (21) 

Here A*",, = 3 is the color factor, and /S^ = (1 — Am'/sy/' the velocity of the final quarks. In 
eq. (|2T|) the leading electroweak radiative corrections have been included by introducing the 
quantities Dajs, a, P = 7, Z, which are defined in terms of the Fermi coupling constant Gp 
and the running QED coupling a{s): 



D 
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ATca'is) 

Dz-y = 



, Dzz 



G'p 



96n {s - M|)2 + (sTz/Mz)' 
M'zis-M'z) 



3V2 {s-M^z? + {sVz/MzY' 



[22] 



At 0(0^), the deviations of the total cross section and the forward-backward asymmetry 
from the tree level values, 5a = a — and 5ApB = ApB — A^p^, are then 



6a =N,(3g [D,,el [(3 - Pl)e,5V^ - fy,S^ 



+Dzz{v! + al) \{3 - fS'MVf + 2(3\6A^ - p^v^S^ 



+Dz-yeeVe 



^2 



^-Ji(vq5V^ + e,5V/) + /?Ja,M^ - ^(e,^f + VqS:^ 



0' 
_£ 

2 • q ' ~H~ ■ q ) I fq^q^^'-q ^^Q^q ' ^<]^q 

Dzjeea,{eq6A^ + aq6V^ + Vq6Ay) + ADzzaeVeiSVfaq + 6A^Vq) 6a 



6A 



FB 



DzyCeaeCqaq + ADzzdeVedqVq 



a"- 



(23a) 
(23b) 



These expressions have to be supplemented by including the standard QCD corrections; the 
formulae for the cross section and the forward-backward asymmetry in the massive case 



can be found in |T^. In the case of the cross section, one can however use the Schwinger 
formulae PDI, which provide a very good approximation to the exact result; this is done by 
performing the following substitution = 7, Z): 



4a, 

3 TT 

4 a, 

3 TT 



TT 



3+^ 
4 

19 



y 



3 

4^ 

7 



22 

10 5 ^ 2^" 



vr 



3 
4 



(24a) 



. (24b) 



A similarly simple yet very accurate substitution also exists for standard QCD corrections 
to the forward-backward asymmetry 

(n")°Kr^(n")°K)°|i + 



a. 



■(3-/3^)^1^ 



(25) 



Finally, on top of the Z resonance these expressions simplify considerably. Besides the fact 
that only the Z exchange contribution has to be taken into account, one can neglect to a 
good approximation the quark masses (except possibly in the h mass matrix; see below) since 
top decays of the Z boson are kinematically forbidden (to achieve a better precision one can 
eventually include the leading mass effects in the Born term as well as the QCD corrections 
in the case of the bottom quark; see [0]). In this case, the 5*^ terms vanish and the deviation 
of the decay width Fg = r(Z — > qq) and the forward-backward asymmetry Aps from their 
tree level values are simply given by 

5Vq _ Vg6Vg^ + ag6A^ 

5Afb_ ^ v,SA^ + agSVf _ ^ VgSVf + a^SA^ 
^FB aqVq vl + al 



3. Results 

We are now in a position to present some numerical examples. In Fig. 2 we show SUSY QCD 
corrections to the hadronic decay width of the Z boson; the solid (dashed) curves are for hh 
(cc) final states. We have set the A parameters in the squark mass matrices (^) to zero, and 
have assumed equal SUSY breaking masses for all squarks, denoted by {friq). Moreover, in 
this figure we have assumed that all parameters entering the squark mass matrices, as well 
as the gluino mass, can be varied independently ("global SUSY" scenario). The four upper 
curves are for negligible mixing between L and R squarks. Even in this case the "D-terms" 
[Dz in eqs.(l)) lead to nonnegligible mass splitting between squarks of different flavour, if 
tan /? 7^ 1. In particular, for tan [3 > 1 (which is favoured by supergravity models 0), m type 
squarks are lighter than d type squarks; as a result, for a given value of {m^ the corrections 
to cc production are larger than those to hb production. 



The uppermost curves in Fig. 2 have been obtained by chosing a very small gluino mass, 
Mg = 3 GeV. A gluino of this mass could have escaped all experimental searches, provided 



squarks are heavier than 100 GeV or so ||23|. Although a careful study showed p4| that a 



GeV gluino does not reduce the slight discrepancy between values of as extracted from low 
energy experiments and those derived from event shape variables measured at ^ M^, 
present measurements cannot exclude its existence, either. It should also be noted that 
squark mass bounds from hadron colliders might be invalidated by the presence of such 
a light gluino. This is because in this scenario, squarks predominantly decay into gluinos, 
which lose a considerable fraction of their energy in QCD radiation prior to their decay, 
thereby leading to a rather soft missing pt spectrum [^. From Fig. 2 we conclude that 



1-loop SUSY QCD corrections to the hadronic width of the Z boson could amount to about 
0.3%, or about 8% of the standard QCD correction. For very light gluinos, 2-loop SUSY 
QCD corrections are also not entirely negligible M]; they amount to about —2% of the stan- 



dard QCD corrections [Q. Altogether SUSY QCD corrections to Fhad therefore amount to 
at most +6% of the standard QCD corrections, for a gluino mass of a few GeV and squark 
masses around 100 GeV. In this scenario the value of as extracted from the measurement 
of Fhad would therefore have to be reduced by about 6%. At the same time, in the presence 
of light gluinos the value of derived from event shape variables has to be increased by 



about 8% 1 24]. The net result is that the present small discrepancy between these two de- 



terminations of as is diminished^. 

If we chose gluino and squark masses above the region excluded by hadron collider searches 
[H the maximal size of the corrections to Z partial widths drops by about a factor of 2, 
as illustrated by the curves for = and Mg = 160 GeV in Fig. 2. Moreover, squark 
mass sphtting due to D-terms becomes less important, so that to good approximation the 
simplified expressions (0) and (pOD can be used. 



Finally, the lowest curve in Fig. 2 demonstrates that squark mixing can have sizable effects 
already for b squarks. The off-diagonal elements of the b squark mass matrix (^) can be 
substantial if tan/3 ^ 1 and n is not too small. Indeed, for the parameters chosen in Fig. 2, 
the lighter b eigenstate would be lighter than 45 GeV, in violation of LEP bounds 0, un- 
less (mq) > 180 GeVQ. In this scenario the correction to the partial width into bb pairs is 
negative. The corrections to the u, d, s, c partial widths for the same set of parameters are 
still positive, however, leading to a very small correction to the total hadronic width of the 
Z boson. In order to test this scenario experimentally one would thus have to measure the 
bb cross section with a precision of a fraction of 1%, which appears to be quite difficult. 

In Fig. 3 we plot the correction to the total cross section for the production of light quarks at 
■\/s= 500 GeV. In this figure we have switched off both squark mixing (by setting Ag = fi = 0) 
and squark mass splitting through D-terms (by setting tan/5=l); however, the previous fig- 
ure showed that results for nonzero and tan/5 7^ 1 are quite similar unless the gluino 

^Strictly speaking the analysis of Ref. is valid only for very heavy squarks; however, we expect it to 
hold also for squark masses around 100 GeV, since squark exchange diagrams contributing to qqgg production 
are not enhanced by large logarithms, unlike diagrams where a gluino pair is produced from a gluon. 



^Note that hadron collider data do not exclude the existence of a single light squark species, if the mass 
of the lightest neutralino exceeds about 15 GeV |p6|| . 



is very light or tan/3 is very large. We see that the corrections reach a maximum at 
rriq ~ OA^/s ~ 200 GeV, almost independently of the value of Mg. If both rriq and Mg 
are much smaller than ^/s the corrections become negative; in the limit of exact SUSY 
{rriq — > rriq, Mg 0) one encounters logarithmic infrared divergencies. For rriq > 200 GeV 
the size of the corrections decreases rapidly. However, even if squarks are not accessible to 
the accelerator we study, i.e. for rriq > v^/2, the correction can be as large as +1%, or about 
one third the standard QCD correction. We also observe that the correction depends less 
sensitively on the gluino mass than on the squark mass; this has also been found in Ref . [|10| . 



In figs. 4a,b we present results for SUSY QCD corrections to tt production at = 500 
GeV, for rrit = 150 GeV. We have fixed the gluino mass to 250 GeV and chosen tan/? = 2. 
The dashed curves are again valid for a "global SUSY" model, with m^^ = m^^ = At and 
jj, = 500 GeV. In contrast, the solid curves are for a supergravity scenario, where scalar 
masses are assumed to be equal to each other, and also equal to the scalar trilinear inter- 
action parameters Aq, at the scale of Grand Unification, Mx — 10^^ GeV. The parameters 
at the weak scale have then be computed by solving a set of coupled renormalization group 
equations [Q; for simplicity we have treated them using the analytical approximations given 



in Ref. |2^. Notice that in this model is no longer a free parameter, but determined by the 
requirement of correct SU(2)xU(l) symmetry breaking, Mz = 91.1 GeV. Moreover, m^^ is 
considerably smaller than m^-^ at the weak scale, due to quantum corrections involving the 
t quark Yukawa coupling. 



The results of Fig. 4 are presented as a function of the mass of the lighter t eigenstate. The 
mass of the heavier eigenstate varies between 390 and 620 GeV in the global SUSY model, 
and between 750 and 1100 GeV in the SUGRA scenario we are considering. Moreover, in 
the former case the t mixing angle is close to 45°, since the diagonal elements of the t mass 
matrix ( |5a| ) are almost equal, while in the latter case the angle is considerably larger than 
45°, so that the light eigenstate is predominantly tn. We see that the correction to the 
total cross section, shown in Fig. 4a, is not very sensitive to the differences between the 
two models we are studying. The reason is that the total cross section is dominated by the 
photon exchange contribution, which does not depend on t mixing, see eqs.(9a) and ( |23a| ). 
The corrections are smaller than for the production of light quarks (with rriq = mi^ ) since 
in case of tt production practically only one squark contributes in the loop, the heavier i 
eigenstate being much more massive. 



In contrast, the forward-backward asymmetry (|23b| ) is sensitive to the Z exchange contribu- 
tion, and hence to t mixing; Fig. 4b shows that for small m^^ even the sign of the correction 
differs for the two models. Unfortunately the absolute value of this correction is always 
less than 0.5%; one would probably need a dedicated "top factory" to achieve this level of 
precision. SUSY QCD corrections to the forward-backward asymmetries of light quarks are 
always well below 0.1%, and can therefore safely be neglected. 



4. Summary and Conclusions 



In this paper we have presented exphcit expressions for the ^qq and Zqq vertices, allowing 
for mixing between the superpartners of left- and right-handed quarks as well as for unequal 
squark masses in the loop. We have found corrections to the total cross section (or, on the 
Z pole, to the hadronic decay width of the Z) to be usually positive, unless both squark 
and gluino masses are much smaller than the centre-of-mass energy a/s- In the limit of 



no squark mixing and degenerate squark masses we reproduce the results of Ref. |IO[.For 
massless quarks, corrections are largest if ~ 0.4-yi, i.e. just above the threshold for open 
squark production, where they can reach +2%; they fall below 1% at niq ~ 0.6-\/i, the exact 
value depending on the gluino mass. If the gluino mass is just a few GeV and squark masses 
are around 100 GeV, which still appears to be allowed experimentally, supersymmetric QCD 
corrections might help to improve the agreement between values of derived from event 
shape variables at ^/s ~ Mz and from the total hadronic decay width of the Z boson. We 
also found that h mixing can be important, and can even flip the sign of the correction. 

The corrections to the total ti production cross section are usually smaller than for the case 
of light squarks, for a given mass of the lightest squark eigenstate of a given flavor. The rea- 
son is that t mixing pushes the mass of the heavier t eigenstate to such large values that its 
contribution is essentially negligible. We also computed corrections to the forward-backward 
asymmetry, and found them to be well below 0.1% for light quarks. In case of t quarks these 
corrections are sensitive to the details of t mixing, unlike the total ti cross section; however, 
even for t quarks the forward-backward asymmetry is changed by less than 0.5%. 

We conclude that, barring the existence of a very light gluino, supersymmetric QCD cor- 
rections to the production of qq pairs in e^e~ annihilation are probably only observable at 
energies above the open squark threshold. They will therefore not be useful as a tool to 
search for supersymmetry; however, after the discovery of a "new physics" signal they might 
allow to confirm its interpretation in terms of supersymmetry. 
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Appendix: Scalar Loop Integrals 

In this Appendix we collect expressions that allow to evaluate the loop functions that appear 



in sec. 2. The scalar one, two and three point functions, Aq, Bq and Co are defined as |18 

(27r/i)"-4 



Ao(mo) 
Bo{s, nil, 1712) 
Co(s,mi,m2,m3) 



171^ 

(27r/i)'^-4 



k"^ — ml + ie ' 



dJ'k 



(27r/i) 



71-4 



ijx? — m\ + ie)[{k — qY — ml + ie] 

dt^k 



(A.l) 



Here n is the space-time dimension and /z the renormalisation scale. 

After integration over the internal momentum k, the function Aq is given by: 



Ao(mo) =ml[l + Ao] 



2 LL^ 

7i5 + log(47r) + log-^, 



(A.2) 



where 7^; is Euler's constant. The function Bq and its derivative with respect to s, -Bq, are 
given by 



„ , , 1 , ^ ^ , ^ mi — mi mi x+ — x_ , X- 

Bo{s, mi,m,) = -(Ai + A2) + 2 + ' log ^ + log — , 



So(s,mi,m2) = - — 



2 + 



mS — m 



-log-i + - 



2 (mf — m\y — s{m\ + 7722) x 



m2 



a:+ — 



log — 



with 



a;± = s — — 777.2 ± y — 2s(77if + 7772) + (^nf — m 



2\2 



(A.3) 



(A.4) 



Note that the x± can be complex. For (7711 — 7772)^ < s < (ttii + 777-2)^, the logarithms 
appearing in eqs.( |A.3|) can be expressed in terms of an arctan of a real argument. When 
writing these equations we have ignored the imaginary parts of Bq and B'^; they are not 
relevant for us, since to next-to-leading order we are only interested in the interference 
between the (real) tree-level and one-loop apmlitudes. 

In this paper we need the three point scalar function Co only for pf = P2 = in this case 
it can be written in integral form as 



Co(s, 777q, 7771, 7772, 7773) 



dy I dx ay + bx + cxy + dy + ex + f 
JO L 



(A.5) 



where 



a = m^ , b = s , c = —s , d = ml — ml — m'^ , e = m\ — m\ , / = 7773 — ie. (A. 6) 



Co can be expressed in terms of a sum of Spence functions Li2(x) = — Jq dt\og{l — xt)/t: 



Co(s,mq,mi,m2,m3) 



5 i=ir- 



Li2 



Li2 



- 1 



,(A.7) 



where we have defined 



Xi 



X2 



2d + e{l-f3,) 1 



2d + e(l - 



l/3± 



-c - e ± y'(c + e)2 - 46(a + + /) 
2b ' 
-d-e± ^{d + eY -Af{a + h + c) 

2(a + 6 + c) 
-d ± Vrf^ - 4a/ 



2a 



(A. 



The and yi± can again be complex. Eq. ([A .71) is only valid above the gg threshold, i.e. for 
s > 4m^. Below the threshold analytical continuation of complex logarithms requires the 
introduction of additional terms; see Ref. [p^ . 
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Figure Captions 



Fig. 1 Vertex (la) and self-energy (lb) corrections to e+e — > qq from supersymmetric QCD. 

Fig. 2 Supersymmetric QCD corrections to the decay width of the Z boson into cc (dashed) 
and bb pairs (sohd). We have assumed a "global SUSY" scenario, where all parameters 
of the squark mass matrices can be varied independently. For c and b squarks the A- 
terms are always negligible; the other parameters are as indicated in the figure. Notice 
that {rriq) is the common SUSY breaking diagonal squark mass, which is also the 
average first generation squark mass, since D-term contributions cancel after summing 
over a complete generation. The lowest curve ends at (m^)=180 GeV since for even 
smaller values, m^^ < 45 GeV. 

Fig. 3 Supersymmetric QCD corrections to the total cross section for the production of light 
qq pairs at an e+e^ collider with y/s=500 GeV. In this figure we have switched off 
squark mixing and chosen tan (3—1, so that the superpartners of all 5 light quarks have 
the same mass rriq. 

Fig. 4 Supersymmetric QCD corrections to the total cross section (a) and forward-backward 
asymmetry (b) for tt pair production at a 500 GeV e+e~ coUider, as a function of the 
mass of the lighter i eigenstate. The dashed curves are for a "global SUSY" model 
with ^ = 500 GeV, while the solid curves are for a supergravity scenario with radiative 
symmetry breaking, where is a derived quantity, as described in the text. In the 
former case we have assumed rrif^ = rrif^ = at the weak scale, while in the latter 
case these relations are vahd only at the Grand Unified scale. 



